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ABSTRACT 



We investigate generation of electromagnetic radiation by gravitational waves interacting with a 
strong magnetic field in the vicinity of a vibrating Schwarzschild black hole. Such an effect may play 
an important role in gamma-ray bursts and supernovae, their afterglows in particular. It may also 
provide an electromagnetic counterpart to gravity waves in many situations of interest, enabling easier 
extraction and verification of gravity wave waveforms from gravity wave detection. We set up the 
Einstein-Maxwell equations for the case of odd parity gravity waves impinging on a static magnetic 
field as a covariant and gauge-invariant system of differential equations which can be integrated as an 
initial value problem, or analysed in the frequency domain. We numerically investigate both of these 
cases. We find that the black hole ringdown process can produce substantial amounts of electromagnetic 
Oh' radiation from a dipolar magnetic field in the vicinity of the photon sphere. 
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C*") ■ 1. Introduction 

J^. i In recent years there has been an enormous effort worldwide to detect gravitational radiation (see, e.g., Barich 

and Weiss (1999); Willke et al. (2002); Ando et al. (2002); Acernese et al. (2002)). It is hoped that within the next 
(T) , few years these detectors will be able to consistently detect and measure the gravity waves (GW) emitted from 
such events as black hole (BH) merger (Buonanno 2002) and exploding and collapsing stars. A pressing problem 
for these detectors is the extraction of the actual waveform from the huge amount of noise invariably generated in 
the detection process. The race is currently on to calculate these waveforms in every conceivable situation in order 
that gravity wave signatures can eventually be statistically extracted from the noise continuously generated by 
these detectors (Flanagan and Hughes 1998a, b; Nicholson and Vecchio 1998), a formidable task. We discuss here 
Q-f a mechanism describing how many of these events could be accompanied by an electromagnetic (EM) counterpart 
q with the same waveform, which could considerably aid in this process. 

Many events will be accompanied by an optical counterpart, such as in supernovae (SN) II and some compact 
binary mergers (Sylvestre 2003), but many in general will not, such as BH-BH merger, and BH ringdown. In any 
case these will only tell us to expect detection, and not the precise form of the waveform to try to extract. What 
would be highly useful for GW detection would be a simultaneous optical detection of the event with the EM 
waveform mirroring that of the GW. This is the situation we discuss here. 

When a plane gravity wave passes through a magnetic field, it vibrates the magnetic field lines, thus creating 
EM radiation with the same frequency as the forcing GW, an effect which has been known for some time (see, e.g., 
Cooperstock (1968); Gerlach (1974); Marklund, Brodin, and Dunsby (2000) and references therein). This would 
provide exactly the mechanism required: virtually all stars have a strong magnetic field threading through and 
surrounding them, and this field becomes immensely strong as the field lines are compressed as the star collapses 
to a BH or neutron star; anything up to 10 14 G - possibly even higher - seems possible in magnetars (Kouveliotou 
et al. 1998). 

It has been proposed that this mechanism may indeed have been observed, being partly responsible for the af- 
terglow observed in some gamma-ray bursts (GRB) and SN events, an argument strengthened by certain anomalous 
GRB and SN light curves (see Mosquera Cuesta (2002) and references therein for a detailed discussion). The basic 
idea is that these events are thought to form a BH or neutron star after the initial explosion (envelope ejection), 
surrounded by a thin plasma which can support a strong magnetic field able to reach supercritical values over a 
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relatively long period of time (compared to the period of the emitted GW). The formation of the compact object 
will release a substantial fraction of its mass as GW which could then be converted into EM radiation as it passes 
through the plasma. Individual models differ considerably in many respects; in particular, additional GW may be 
produced from a stressed accretion disk powered by the spin of the BH (van Puttcn 2001). 

Studies of generation of EM radiation by GW in astrophysical situations so far have provided order of mag- 
nitude estimates (Mosquera Cuesta 2002) and some of the extra complexities involved when a thin plasma is 
present (Macedo and Nelson 1983; Daniel and Tajima 1997; Brodin and Marklund 1999; Marklund, Brodin, and 
Dunsby 2000; Servin, Brodin, and Marklund 2001). In particular, a thin plasma can increase the frequency of the 
electromagnetic radiation, whose origin is from a plane gravity wave passing through a uniform, static, magnetic 
field, thus strengthening the observational potential of the EM-GW interaction still further (Brodin, Marklund, 
and Servin 2001; Servin and Brodin 2003). While these investigations have given a good indication of the physical 
processes we may expect, the effect has not yet been studied in a strong gravitational field, the most promising 
place we may expect such an interaction as likely to happen. 

Our aim here, then, is to study the induced EM field from the interaction of GW emitted during BH ringdown, 
the settling down of a BH after an initial perturbation, with a strong magnetic field which surrounds the BH. Shortly 
after a BH is disturbed by any kind of small perturbation, it radiates its curvature deformations as GW with certain 
characteristic frequencies which are independent of the initial perturbation, and dependent only on its mass (in 
the case of a Schwarzschild BH, which we consider here). These complex quasi-normal frequencies form solutions 
known as quasi-normal modes which govern the BH ringdown process (Nollert 1999; Kokkotas and Schmidt 1999). 
As the ringdown process is thought to be independent of the initial perturbation, we may expect that studying 
this particular situation will help give understanding in more complex situations, such as the late stages of BH-BH 
merger (Buonanno 2002). Indeed, while it would seem logical that perturbations of Schwarzschild would give little 
information about something as non-linear as colliding black holes, it turns out that perturbation theory gives 
surprising accuracy in many cases of interest (Price and Pullin 1994). 

The frequency of this generated EM radiation will be very low, generally less than about 100kHz, and would 
be typically absorbed by the interstellar medium. This is where the photon frequency conversion (Marklund, 
Brodin, and Dunsby 2000; Brodin, Marklund, and Servin 2001; Servin and Brodin 2003) could come into play, 
overcoming this by increasing the frequency to detectable levels. An important extension of this work, therefore, 
will be to include a plasma into this situation; we leave this to later, and concentrate here on setting up a suitable 
formalism for the inclusion of a plasma, while investigating the pure curvature effects of the BH, which turn 
out to be quite large. We find that the amplification of the EM field is much stronger than is the case of plane 
GW (Marklund, Brodin, and Dunsby 2000), where the amplification grows linearly with interaction distance. Here, 
we find substantial growth in the vicinity of the horizon and photon sphere. 

1.1. An overview 

Electromagnetic waves around a Schwarzschild BH generated by gravitational waves interacting with a strong 
static magnetic field are governed by the Einstein-Maxwell equations which are of the form 

EM field around a BH = (GW x Strong static magnetic field) . (1-1) 

The homogeneous solution to these equations, where the induced currents are zero, will be governed by the well 
known Regge- Wheeler (RW) equation for EM perturbations around a BH (Price 1972), while the GW terms are 
governed by the Regge- Wheeler equation for gravitational perturbations of a BH (Regge and Wheeler 1957). The 
RW equation describes how the fields of different spin are lensed and scattered around the hole. A description 
of the more general situation including the rhs of Eq. (1-1) is rather less trivial than Eq. (1-1) may suggest, for 
several reasons: for example, the GW terms on the rhs need some manipulation to convert them to the familiar 
RW variable, and the existence of the magnetic field on both sides of the equation means that gauge problems 
are paramount, and considerable work must be done to cast the equations into a manifestly gauge-invariant form. 
While it may be expedient to use the Newman-Penrose formalism (Chandrasckhar 1983) for this problem as all 
variables in the perturbed spacetime travel on null cones, an important extension of this work will be to include 
plasma effects, to model a more realistic astrophysical environment, for which the Newman-Penrose formalism is not 
so adept. In addition, the electric and magnetic fields require a timelike vector field for their definition. For these 
reasons we will use the covariant and gauge-invariant perturbation method introduced in Clarkson and Barrett 
(2003), which is ideally suited to this particular problem involving spherical symmetry, and has the advantage that 
is is well adapted to fluids for later use when investigating plasma effects. 

An important issue in rclativistic perturbation theory is the mapping, or gauge choice, one makes between 
the background and perturbed model; many perturbation approaches are not invariant with respect to this gauge 
choice. Metric based approaches to perturbation theory suffer from this gauge freedom, whereby spurious gauge 
modes exist and must be identified. While these spurious gauge modes may be eliminated in analytical treatments, 
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when the equations are integrated numerically these modes have a tendency to grow very fast without bound, in 
so-called gauge pathologies (Alcubierre 1998). Furthermore, the tractability of the problem often depends upon 
judicious gauge choice; hardly an ideal situation (Ruoff, Stavridis, and Kokkotas 2001). 

The covariant '1+1+2' approach we utilise relies on the introduction of a partial frame which form the differ- 
ential operators of the spacetime, and allow all objects to be split into invariantly defined physical or geometric 
objects. Covariant perturbation techniques initiated in Ellis and Bruni (1989) are employed to write the equations 
in a fully gauge-invariant form which can then be solved with the use of appropriate harmonic functions which 
remove the tensorial nature of the equations. To aid in the solution we consider it formally as a second-order 
perturbation problem, and introduce 'interaction variables' for quadratic quantities. This then allows us to write 
the equations for the induced EM radiation as a system of gauge-invariant, covariant, first order ordinary differ- 
ential equations in the relevant variables, while we can easily convert to covariant wave equations for clarity and 
integration as an initial value problem when desired. 

The paper is organised as follows. The covariant formalism we use is reviewed in Sec. 2. In Sec 3 we derive 
the coupled perturbation equations which govern the interaction. These equations are integrated numerically in 
Sec. 4, and the implications for the emitted radiation discussed. We briefly conclude in Sec. 5. An Appendix gives 
some key formulae relating to the spherical harmonics we use. 

Sections 2 and 3 contain most of the technical material, the crucial results being Eqs. (3-41) - (3-43) and 
(3-45); some readers may wish to skip to Section 4 where specific astrophysical situations are discussed. 



2. The 1+1+2 covariant approach 

Covariant methods in General Relativity (GR) are formulated in a very different way from coordinate metric 
based approaches: in the latter, Einstein's field equations (EFE) are second order partial differential equations 
in the components of the metric; in the former, a physical (partial) frame is chosen and the Ricci and Bianchi 
identities are irrcducibly split with respect to this frame, resulting in a system of first-order differential equations. 
Supplemented by the crucial commutation relations between the frame vectors, this system of equations becomes 
equivalent to the EFE, but always deals with invariantly defined physical or geometric quantities. 

The 1+1+2 covariant sheet approach relics on the introduction of two frame vectors: the first being a timelikc 
vector field u a : u a u a = —1, representing the congruence on which observers sit; the second is a spacelike vector 
field n a : n a n a = 1, n a u a = 0, which can be chosen along a preferred direction of the spacetime. These two vector 
fields define a projection tensor 

N b = h b -n a n b = g b + u a u b -n a n b , (2-1) 

which projects vectors orthogonal to n a and u a : n a N a b = = u a N ab , onto 2-surfaces (N a a = 2) which we refer to 
as the 'sheet'; h a b = g a b + u a u b is the tensor which projects orthogonal to u a , into the observers' rest space. Using 
h a b, any 4-vector may be split into a (1+3 scalar) part parallel to u a and a (1+3-vector) part orthogonal to u a . 
Any second rank tensor may be covariantly and irrcducibly split into scalar, 3-vector and projected, symmetric, 
trace-free (PSTF) 3-tensor parts, which requires the alternating tensor e a bc — u d r)d a bc', these are the key quantities 
in the 1+3 covariant approach (Ellis and van Elst 1998). Crucially, the covariant derivative of u a may be split in 
the standard manner, the irreducible parts (the acceleration ii a , the expansion, 6, the shear, a a b, and the rotation, 
u> a ) forming some of the key variables of the 1+3 approach - we refer to Ellis and van Elst (1998) for further details. 

Now, using N a b, any 3-vector tp a can now be irreducibly split into a scalar, ^, which is the part of the vector 
parallel to n a , and a 2-vector, lying in the sheet orthogonal to n a : 

iP a = ^ n a + #°, (2-2a) 

where 

* = ip a n a , and 5< a = N ab ^ h = V, (2-2b) 

where we use a bar over an index to denote projection with N a b. Similarly, any PSTF tensor, ip a b, can now be 
split into scalar, vector and 2-tensor (which are PSTF with respect to n a , and is therefore transverse- traceless) 
parts: 

4>ab = V><a&> = * (n a n b - \N ab ) + 2* (a n b) + * afc , (2-3a) 

where 

* = n a n b ip ab = -N ab ^ ab , (2-3b) 
* a = N a b n c Ac = * a , (2-3c) 
* ab = ^ {ab} ee (N (a c N b) d - ±N ab N cd ) Ad- (2-3d) 
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We use curly brackets to denote the transverse-traceless (TT) part of a tensor. We also define the alternating 
Levi-Civita 2-tensor (area element) 

£ab = e abc n c = u d i] dabc n c , (2-4) 

so that e a bn b = = £( a b)- 

With these definitions, then, we may split any object into scalars, 2-vectors in the sheet, and transverse- 
traceless 2-tensors, also defined in the sheet. These three types of objects are the only objects which need to be 
solved for, after a complete splitting. Hereafter, wc will assume such a split has been made, and 'vector' will 
generally refer to a vector projected orthogonal to u a and n a , and 'tensor' will generally mean transverse-traceless 
tensor, defined by Eq. (2-3). 

We split the familiar 1+3 variables in this manner; in particular, the electric and magnetic fields are irreducibly 

split 

E a = Sn a + £ a , (2-5a) 
B a = @n a +3S a , (2-5b) 

while the kinematical and gravitational variables become, using Eqs. (2-2) and (2-3) 

ii a = An a + A a , (2-6a) 

u a = nn a + n a , (2-6b) 

a ab = S (n a n b - jN ab ) + 2S (a n b) + S a6 , (2-6c) 

E ab = £ (n a n b - \N ab ) + 2£ {a n b) + £ ab , (2-6d) 

H ab = U (n a n b - \N ab ) + 2H (a n b) + U ab . (2-6e) 

For example, £ — n a n b E ab is the tidal force along n a , £ a = N a b n c E bc is a 'drift' vector, while £ ab = E^ ab j is the 
TT part of the electric Weyl curvature in the sheet orthogonal to n a . 

There are two new derivatives of interest, which n a defines, for any object ■0... " : 

^ a ... b c - d = n^ a ... b c - d , (2-7a) 
6 e ^ a ... b c - d = NjNj ■ ■■N b g N h ° ■ ■ ■N i d D j ^ f ... g h - i , (2-7b) 

where D a is the spatial derivative defined by h b (see, e.g., Ellis and van Elst (1998)) . The hat-derivative is the 
derivative along the vector field n° in the surfaces orthogonal to u a . It is important to note, however, that these 
derivatives do not commute; commutation relations for scalars are given in Clarkson and Barrett (2003). This is a 
vital aspect of the formalism. 

With these definitions we may now decompose the spatial projection of the covariant derivative of n a orthogonal 
to u a : 

D n 6 = n a a b + \(j>N ab + £e ab + ( ab , (2-8) 

where 

a a = n c D c n a = n a , (2-9a) 
^ = S a n a , (2-9b) 
£ ee \e ab 5 a n b , (2-9c) 

Cab = 5 {a n b} . (2-9d) 

We may interpret these as follows: travelling along n a , <p represents the sheet expansion, Q ab is the shear of n a 
(distortion of the sheet), and a a its acceleration, while £ represents a 'twisting' of the sheet - the rotation of n a . 
The other derivative of n a is its change along u a , 

h a = Au a + a a where a a = n a and A = n a ii a . (2-10) 

The new variables a a , <p, £, ( ab and a a are fundamental objects in the spacetime, and their dynamics gives us 
information about the spacetime geometry. They are treated on the same footing as the kinematical variables of 
u a in the 1+3 approach (which also appear here). 

The 1+1+2 split of the Ricci identities for u a and n a , and the Bianchi identities, provide a complete set of 
first-order differential equations for these variables, and were discussed in Clarkson and Barrett (2003) for the case 
of a gravitationally perturbed BH. We shall not require a further generalisation of these equations here. 
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Obviously, wc shall require Maxwell's equations (ME), which may be irreducibly split using these definitions: 

g + 5 a S a = -<K + S a a a + 2Q.3S + m a SS a + /j, p e , (2-lla) 

J> + 5 a S8 a = -4>S8 + SS a a a - 2Qg - 2fl a g a , (2-llb) 

g - e ab S a ^ b = 7£,& + S a a a - (|6» - S) g + £ a <f a + e ab {A a SS b + n a g b ) - /j, J, (2-llc) 

38 + e ab S a g b = -2£g + 3§ a a a - (§0 - E) & + Z a 33 a - e ab (A a g b - n a 33 b ) , (2-lld) 

g- a + s ab ( J" - & b d) = z& a - (|0 + 4) e a6 ^ fc - (fe + ±e) ^ - r» £ab( r b 

+g (-a a + S a + £a6 6 ) + 38e ab {A b - a b ) + Y> ab S b - e ab ( bc 3§ c - fi Ja, (2-llc) 
3B- a - e ab (i b - S b g^j = -^g a + (10 + A) e ab S b - (p + |S) - r!e ab ^ fc 

+^ (-a„ + S a + e ah r! b ) - Se ab (A b - a b ) + S ab ^ b + e ab ( bc g c . (2-llf) 

Here, MKS units are used (/xo), p e is the charge density, and the current density j a has been split into its 1 + 1 + 2 
parts, J and J a . The first two equations arise from the constraint ME, while the rest are the evolution ME. In 
flat space in the absence of currents and charges the rhs of these equations vanish (for a static 'natural' choice of 
frame). Thus, gravity modifies ME in the form of generalised currents. Note how the rotation terms £, fl and £l a 
flip the parities of the EM fields (e ab is a parity operator - see the Appendix). 



3. Electromagnetic radiation around a vibrating black hole 

There are two ways to proceed in solving this problem, depending on how one views Maxwell's and Einstein's 
equations. If we view Maxwell's equations as being essentially separate from the field equations, deciding on the 
fly whether to include the gravitational effects of the EM field, then this particular situation may be considered as 
having the EM field as a test field on a vibrating BH background. An alternative viewpoint is to consider Maxwell's 
and Einstein's equations as a coupled system of equations (with /is ~ \B 2 , etc.), with decoupling occurring only 
when one can legitimately set terms 0(B 2 ) to zero, clearly more intuitive in a perturbation approach. They are 
mathematically equivalent in vacuo only due to the linearity of ME, a feature not present in plasmas in general. 1 
By using the second interpretation things will automatically be easier when complicated non-linear plasma effects 
are included at a later date. We will therefore treat this as a perturbation problem at second- order in a two 
parameter 'expansion' in two 'smallness' parameters egg representing the magnitude of the static magnetic field, 
and e g representing the amplitude of the GW (these are labels for the two types of first-order perturbations as much 
as anything else; we need both because we keep terms 0(eage g ), but neglect terms 0(e 2 ) and O(e^)) - see Bruni 
et al. (1997); Bruni and Sonego (1999). Consequently, this expansion allows us to set up the equations as a system 
of linear first-order differential equations at second-order in the perturbation, which at the same time serves to 
illustrate a new technique for covariant and gauge-invariant non-linear perturbation theory. Once a gauge-invariant 
formalism has been set up to study this interaction it should then be relatively easy to include complicated plasma 
effects and so on. 

We will divide up the perturbation 'background' spacetimes and denote them as follows: 

- B = Exact Schwarzschild, O(e ); 

- T\ = Exact Schwarzschild perturbed by a pure static magnetic field, neglecting the energy density of the 
field in comparison to the curvature of the BH: <D{egg); 

- = Schwarzschild with gravitational perturbations 0(e g ), as given in Clarkson and Barrett (2003); Regge 
and Wheeler (1957); 

- S = T\ + T2 allowing for interaction terms in Maxwell's equations: the induced EM fields will be C(e^e 9 ); 
this is the situation of interest. 



^^We also consider the second interpretation more appropriate in a covariant approach because the frame derivatives in ME auto- 
matically couple the curvature of the spacetime to the EM field through the commutation relations [for example, D^ a D b ^V c brings in 
both Ricci and Weyl curvature in general, and we use this commutator after defining gauge-invariant variables - see Eq. (3-36) below]. 
On the other hand, the coupling to gravity which occurs explicitly in ME in any curved spacetime is from purely kinematical quantities 
(frame motion), which then couple to gravity through the field equations; this is why gravity can be explicitly 'decoupled' from the 
vacuum ME. The commutation relations thus confuse this issue, because through these relations curvature directly induces the EM 
field. 
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We will generally refer to terms of order 0(e@) and 0(e g ) appearing in T as 'first-order', while those variables 
of order C(e^e g ) in S as 'second-order' variables. If one prefers to view ME as a test field (where the only 
'perturbation' is in the EM field) , then these backgrounds may instead be thought of as useful labels for each type 
of field present. 



3.1. The background fields 

We now review each of the backgrounds. 

B: the exact SCHWARZSCHILD solution. For a family of static observers, in the background we have only the 
zeroth-order scalars: £, the radial tidal force; A, the acceleration a static observer must apply radially outwards (to 
prevent infall); and <fi, the spatial expansion of the radial vector n a . These are determined by the radial propagation 
equations 2 : 



-¥ 2 



£, (3-la) 

£ = -§<££; (3-lb) 

together with 

£ + A(f> = 0. (3-lc) 
Defining the affinc parameter by " = d/dp, and another radial parameter r by 



r = 



\<fr, (3-2) 



the parametric solution to these equations, giving a complete description of the BH, are given by 



£ = 


2m 




<t> = 




2m 






r 


A = 


m ( 
~, (I" 


2m 
r 



-1/2 



where 



r \ L 2m 



(3-3a) 
(3-3b) 

(3-3c) 



p = 2mcosh-^-j + ^l--, (3-4) 

relates the affine parameter p associated with the radial vector n a with the usual Schwarzschild coordinate r. 
In T and S we keep all powers of these variables. 

T\: the static magnetic field. The following equations govern the static {38 = 08 a — 0) magnetic field: 

J = -5 a 3g a - 4>B8, (3-5a) 
& a = 8 a &-{\<l> + A)& a , (3-5b) 
= e ab S a ,^ b . (3-5c) 

The last equation tells us that the field is purely of even parity. In general the solution to these equations when 
harmonically decomposed can only be written as a complicated combination of hypergeometric functions (which is 
partly why the perturbation method we utilise below is effective). 
In S we neglect all products of the magnetic field with itself. 

The solution for a dipole field is of particular importance: when split into spherical harmonics (see Clarkson 
and Barrett (2003) and the Appendix), the 1=1 equations have two solutions; one which is uniform at infinity, 
characterised by SS = and one which falls off like 1/r 3 at infinity, which is the true dipole. The solution for the 



2 We refer to equations in involving the radial hat-derivative as propagation equations, and those involving the temporal dot-derivative 
as evolution equations; equations involving neither of these may be thought of as constraints, though this depends on how one chooses 
to integrate the equations - see below. 
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latter part is, in terms of r: 



8m 3 
8m 3 



In 1 



2m 
r 



2m_y 1/2 



2m 
r 

1 - 



2m^ 



2m^ 



+ 



2m 



m\ 
r / 



(3-6a) 
(3-6b) 



where is the magnitude of SS^r 3 as r — > oo. 



•F 2 : THE gravity wave perturbation. As shown in Clarkson and Barrett (2003), these perturbations arc 
governed completely, in the A a = S a (p = 5 a A = frame, by the tensorial form of the Regge- Wheeler equation 
(Regge and Wheeler 1957; Clarkson and Barrett 2003) 



-W ab + W ab + AW ab - <j> 2 W ab + 8 2 W ab = 0, 



(3-7) 



where the Regge-Wheeler tensor W a b is a gauge- and frame-invariant TT tensor, defined as (Clarkson and Barrett 
2003) 



W ab 



2 Cab - \r 2 £- x 5 {a Xb h 



(3-8) 



and X a = 5 a £ is the gauge- invariant variable (Stewart and Walker 1974) describing the angular fluctuation in the 
radial tidal force. This tensor contains in compact form the curved space generalisation of the two flat space GW 
polarizations h + and h x 3 . 

Every other object in Ti is determined by linear combinations of {Wj, Wj, Wj Wt}, once appropriate har- 
monics are used (see the Appendix and Clarkson and Barrett (2003)). While Eq. (3-7) governs GW of both parities, 
for simplicity we shall only consider the case here where the GW are of odd parity. For purely odd perturbations 
the gravitational field is governed by W a b, and the other GW variables that we shall require are related to this by 
the covariant gauge-invariant equations (Clarkson and Barrett 2003) 



Cab = 


2 

-T-oWab, 


(3-9a) 


S a fc = 


r- 2 W ab + ^W ab , 
(pr z 


(3-9b) 


a £a 


^s b w ab , 


(3-9c) 


£ab'H b = 


-^ b w ab , 

<pr A 


(3-9d) 


H = 


--^e ab 5 a S c W b c . 
(hr z 


(3-9e) 



Although they may be given in a similar fashion, we shall not require £ a b,Hab- (All other 1+1+2 variables arc 
zero.) While we would not normally require parts of the Weyl tensor to solve ME, we will need them here as they 
arise when generating propagation equations for the gauge-invariant part of the magnetic field in S, through the 
commutation relations. We shall find that turning ME into a gauge-invariant system at second order explicitly 
introduces Weyl curvature into the problem. 
In S we neglect all products of these quantities. 



Because the background B is spherically symmetric, the solutions of both parts of T may be expanded in 



3 In Minkowski space, in the transverse-traceless gauge, gravitational waves are described by the TT tensor, 

" I h x -h+ J 

(see Box 35.1 in Misner, Thorne, and Wheeler (1973)). In cartesian coordinates for a plane wave travelling along the z-axis, £ MJ/ = 
\dzh^, . Thus, far from a localised source of radiation, the scale-invariant part of the GW is W a b — rQ a b which is the power of the 
GW; hence (W a tW ab ) (where the angle brackets denote the average over several wavelengths) is proportional to the energy carried 
from the source by the GW (Misner, Thorne, and Wheeler 1973). 
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spherical harmonics. This implies that we can write 



SB = J2 ^ ' ( 3 " 10a ) 

= E ^ - (3-10C) 



,=2 



where the g and 23 subscripts serve to remind which harmonic indices we are summing over in each case, a 
distinction required in the next Section. Then the harmonic components of the magnetic field and RW tensor obey 
the constraint equations, where L = I (£ + 1), 



g2£V> = _ L ^ r -2^%) ; (3 _ lla) 

= (1 - L a ) r^df , (3-llb) 

e ab S a ^ } = (3-llc) 

and 

5 2 W% = (4> 2 -3£-L g r-2)W a ^ (3-12a) 

6 a 5 b W^ = 0. (3-12b) 

Assuming separable solutions implies the usual spherical harmonics for the angular parts sec the Appendix. We 
can write the variables in this way because only one parity is present for each field. 



3.2. The interaction terms in Maxwell's equations 

Here we will introduce a set of auxiliary variables, all of order 0{egge g ), which allow us to convert ME into 
a linear (in differential order) system of gauge-invariant ordinary differential equations (gauge-invariant because 
they vanish at all perturbative orders lower than this (Bruni et al. 1997; Bruni and Sonego 1999)). We refer to 
these as the interaction variables 

A quick glance at the rhs of ME, Eqs. (2-llb) - (2-1 If) , reveals that we're dealing with products of tensorial 
spherical harmonics, which are not particularly pleasant. Instead of explicitly using tensor spherical harmonics in 
the GWxB products in ME, we shall absorb them into the following interaction variables, which makes the resulting 
equations considerably neater. There is no extra work involved here, although it may not appear that way; we 
would otherwise still require the key equations (3-15) and (3-24) given below. The latter in particular are crucial 
relations among all the coupled tensor/ vector/scalar spherical harmonics which appear (these are the products 
given by Eqs (3-14), (3-22d), although we have absorbed the magnetic field strength and the GW amplitude). 
There is another reason for defining the variables in the manner we do: while variables such as a a appear in ME, 
our solution in T2 only gives us a a ', we circumvent this problem by absorbing the time derivatives into our new 
variables below. 

With these considerations in mind, we define the four interaction variables 

(%.<y / Ve*tg> < f »y (4»4i>\ ,„ , oN 

Xa = [Xla i X2a i X 3 a i X 4 a ) O 13 ) 

as follows: 

r^^-'C^,, (3-Mb) 

il^ = {^T 1 W^5^, (3-14c) 
X ( r ) = (^ 2 )- 1 C^ ! (3-14d) 
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for each l@ l g interaction. We use a bold font as a matrix shorthand for the '4-vector' these variables form. 
These variables obey the propagation equations 



J e » e g> 
Xa 



L Xa 



(3-15) 



with the interaction matrices given by, for each lag and 



jltse,tg) 



( 



V 



■A) 



L.sgr- 2 




1 1 
• + 2A) 

-2(f) 



(3-16) 



where 

A (<y = -w 2 + 3£ + L g r-' A . (3-17) 

We have introduced the time harmonics of Clarkson and Barrett (2003) into these equations for notational sim- 
plicity; factors of iw just represent time derivatives, d/dr - we will discuss the significance of these later. For now 
note that 

2w> - 1/2 
r 



-Aw 



all 



2a 

(f>r ' 



(3-18) 



arising from the commutation relation between the dot- and hat-derivatives. Here, a is a constant, which we will 
discuss below. 

In order to simplify our presentation, we will define a set of auxiliary interaction variables as follows (some of 
which may be a little surprising, but they are all required). First, define 



/ (%) 
(V I V-, 

\r _ la 

y a — U a ) 



5 b W„ 



Vg) 



V, 



2 a 



5 b W 



(3-19) 



a 6 



where we use a bold font to denote the '2-vector' matrix. Similarly we define 



^la 

A, 



2 a 



(3-20) 



For simplicity of presentation, we introduce the shorthand notation 'o' which takes two '2-vectors' to form a 
'4-vector' as 



V o\ = {V 1 \ 1 ,V 2 \i,V 1 \ 2 ,V 2 \2) ■ 
We use these to define the following '4-vector' variables as follows 



(3-21) 







(3-22a) 


*a 


oh' 1 V a oS C X c , 


(3-22b) 


M a 




(3-22c) 


■(t&lg) 
Ja = C 




(3-22d) 



where, for example, K = (K\, K 2 , K3, K4), gives the shorthand for four of these sixteen new variables. These 
variables are all 0{eage g ). They are all constructed to obey the same propagation equation as x, viz: 

(3-23a) 
(3-23b) 
(3-23c) 
(3-23d) 

We have defined all these variables as the time integral of combinations of the RW tensor and the static magnetic 
field. This is because for some of the GW variables appearing in ME it's their time derivatives which are related 





S ,tg) 




*a 






- (% 
Ma 






Ja 
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to the RW tensor [see, e.g., Eqs. (3-9b) and (3-9c)]. Denning auxiliary interaction variables this way which satisfy 
the propagation equations (3-23) removes this problem, and absorbs it into the initial (or boundary) conditions. 

By taking various <5-derivatives of these variables and using the appropriate commutation relations (see Clark- 
son and Barrett (2003)), together with Eq. (3-12b), we can show that they all obey the following constraints, which 
are crucial identities for consistency of the resulting equations later, and allow us to relate all the interaction terms 
to x when we split ME into spherical harmonics. 

= JT" + + 2re ab S b ^ + (L g - 2) - ^ ' , (3.24a) 

= HP*"* + e ab 6^ eg) - (L g 2) e ab ^ X \^ (3-24b) 

= [rS 2 + (L g L m ) r" 1 ] ^ + e ah b^\^ - Te^M^, (3-24c) 

= Lag [S 2 + (L a - L g ) r- 2 ] - v [S 2 + (L m + L g ) r~ 2 ] e ab 5 b ^ k) , (3-24d) 

= (r 2 5 2 + L^) S^T^ - (L g - 2) L^xT^ , (3-24e) 

= (r 2 5 2 + L g ) S^T® 2L g 5 a MT l9) . (3-24f) 

These twenty-four constraint equations propagate consistently. 

For each tag <-> £ g interaction, the system of equations describing the gravitational wave - magnetic field 
interaction are given above. Not all these variables appear explicitly in ME, but they couple to them through the 
system of propagation equations (3-23) and constraints (3-24). We now discuss how these enter ME. Consider, for 
example, the term ( ab ^ h which appears in the evolution equation for S ai Eq. (2-lle). We can relate this to the 
interaction variables above as follows: using Eq. (3-9a) we have 

= foe) < 3 - 25 > 

Similarly for the other products. We therefore use the following abbreviations: 

Xa =^X^\ M a =^M^\ J =E^ (3-26) 

tst,tg ias,e 3 taJg e&,e 9 

while for \I/ a we define: 

*«=Z^(*la >*2a ,V*3a . V*4« J- (3-27) 

tm,e g 

(The definition for \& is slightly different because it is defined having a S 2 ,^'^ ~ L^Si^ term in it.) These 
definitions prevent summations appearing explicitly later. 

3.3. The gauge- invariant form of Maxwell's equations 

Neglecting terms O (egg x even parity gravity waves) and O(e^) (strictly speaking we are only neglecting terms 
O(e^g), see the end of this section), and choosing the frame in T 2 such that A a — 5 a <p = 5 a A = we find that ME 
become, 

i + b a S a + <j>£ = 0, (3-28a) 
J 1 + 5 a 3S a + 4>@ = 0, (3-28b) 
i - e ab 5 a ^ b = 0, (3-28c) 
38 + e ab 5 a <g b = 0, (3-28d) 

and 

4 + e ab (@ b - 8 b <%) + (±</> + A) e ab % h = -2s ab x\, (3-28e) 

A - e ab (i b - 5 b S) - (±<j> + A) e ab S b = 4* 3 a + 4>Xia + 2 X 2a- (3-28f) 

The terms on the left are those which govern an EM field around a BH; those on the right are the interaction terms. 
Note that these equations are a mixture of first and second order quantities, and are thus not gauge-invariant, and 
therefore not integrable. 
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In order to convert ME into gauge-invariant form, it is not enough to define the interaction variables above; 
we must also do something with the magnetic field: in S the magnetic field appearing in ME has a contribution 
from the static background field in T\ which we must somehow subtract off. The standard route to do this is as a 
series expansion, but this does not work here. If we imagine that B a is written as a power series, 

B a = egg {Bt + e g B a 2 +■■■) (3-29) 

where B\ satisfies the T\ equations, (3-5c), then one would imagine B^ would cancel out of the S ME when B a 
appears alone leaving just £? 2 ; when it appears multiplying an Ti term it is only B\ which contributes. However, 
this is not the case. It is possible to show from the commutation relations for the hat and dot derivatives acting on 
B a that this leads to an inconsistency, implying that the interaction terms must be zero. Consider, for example, 
the scalar part of the magnetic field: 

38 = e®38 1 +e g e®38i + 0(e 2 g ,e%) (3-30) 

where 3§i satisfies 3§\ = and 3S\ = F where F = 0, representing the background solution T\. Now, using the 
commutation relation given by Eq. (30) in Clarkson and Barrett (2003), 

3$ = E g Eog 3$ 2 = e g Eog ( J> 2 - A# 2 ) , (3-31) 

by using the commutation relation after substituting from Eq (3-30), and neglecting terms 0(e g ). Alternatively, 

& = &-A38- 2aJ a 38 = e g e m (j 2 - A# 2 ) - 2e sg aJ a 38 1 , (3-32) 

where we applied the commutator before using the expansion given by Eq. (3-30). This is clearly a contradiction if 
a a 8 a 3§i ^ 0, which is the case here. [The correct form of calculating this equation results in Eq. (3-35).] In fact, 
this problem usually arises when using covariant (partial-)frame methods for second order perturbation theory 
In contrast to metric-based approaches, the solutions for perturbed derivative operators are never sought, so they 
must always operate on quantities of the same perturbative order. We must therefore define some gauge-invariant 
variables for the magnetic field. 



3.3.1. Gauge-invariant variables for the magnetic field 

We define the variables 

= 38, and a = 38- a = 38 a + (a b 38 b ) n a , (3-33) 

which are gauge-invariant in S, as they vanish in T (Bruni et al. 1997; Bruni and Sonego 1999). To convert ME 
into a gauge-invariant system of equations, we must somehow replace every occurrence of 33 with 0, and 33 a with 
a . Note first that 

= -e ab 5 a S\ (3-34) 

immediately from Eq. (3-28d). Meanwhile, the commutation relation between hat- and dot-derivatives (see Eq. 
(30) in Clarkson and Barrett (2003)), when applied to 3$ results in the propagation equation 

0=-{4> + A)0- 5 a a + AS a ^ + cp6 a xt + 2^X2, (3-35) 

where we have used Eq. (3-28b), and the appropriate commutation relation for dot-5 derivatives on vectors. How- 
ever, this equation also arises from propagating (3-34) using ME, as it should. Hence, because Eq (3-34) is a 
consistent constraint, the propagation equation for [3 is redundant. This implies that Eq. (3-34) can replace 
Eqs. (3-28b) and (3-28d). To find a propagation equation for (3 a we must propagate Eq. (3-33) using the appropri- 
ate commutation relation for vectors, giving 

Pa = £ab£ b - {\<j> + 2A) a + 6 a - 2 X i a - <fea - ^Xia - 2r~ 2 * ' la - 4r~ 2 M la - 2r- 2 J la , (3-36) 

which replaces Eq. (3-28e). It is this equation which brings Weyl curvature into ME through the commutation 
relations. A key remaining evolution equation comes from calculating § using Eq. (3-28c): 

S = e ab S a b + 4>e ab 8 a X \ + 2e ab 5 a x b 2 , (3-37) 

which propagates consistently. We will use this evolution equation, which is just the gauge-invariant form of Eq (3- 
28c), to replace Eq. (3-28a). Therefore, ME are now just the two vector propagation equations (3-36) and (3-28f), 
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together with the two scalar non-propagation equations (3-34) an d (3-37). The last two serve as definitions for (i 
and S after time harmonics are used; these then become constraints. 

Note how converting the gauge-dependent form of ME, Eqs. (3-28a) - (3-28f), which contain a mixture of first 
and second perturbation orders, into a gauge-invariant second order system has introduced many more interaction 
terms into the equations, terms arising purely from the Ricci identities. These terms are essentially hidden in the 
frame derivatives (dot, hat and 5) when acting on 3§ in Eqs. (3-28a) - (3-28f), and in 38 itself, illustrating the 
importance of using a full set of gauge-invariant variables. 

Although equations we have derived are gauge-invariant to order 0(cgge g ), they are actually valid up to 
0(e 2 ag€ g ), C(e^), which can be easily seen as follows. If we include terms O (e^) (i.e., the energy density 
and anisotropic pressure of the static magnetic field) in the gravity sector, then changes to W a b (and other GW 
variables) is O (e%), making the change to the GW-B variables O (egg[e g + £%])■ However, the equations are not 
gauge-invariant at this order, because the variables are non-zero at lower perturbative order (i.e., O (eage g )); see, 
e.g., Bruni et al. (1997); Bruni and Soncgo (1999). 

The gauge-invariant form of the equations now shows exactly the terms and couplings involved in generating 
the EM field. Consider, for example, the covariant wave equation for S : 

-S + S + 5 2 S + (20 + A)i- (±</> 2 + 2£) S = -20e ab <5 a x ? - 4e ab 5 a X b 2 - 4£ o6 J°tf §. (3-38) 

The lhs of this equation is just the contribution from the BH geometry, and can be related simply by a change of 
variables to the usual Regge- Wheeler equation for an electromagnetic field around a BH [compare with Eq. (3-7) - 
see also Eq. (3-41)]. The rhs, on the other hand is the source from the interaction terms, and has contributions from 
the time integral and angular derivative of the dot-product between the transverse traceless RW shearing tensor 
and the angular (sheet) part of the magnetic field, and the 2-divergence of the RW tensor times the magnitude of 
the radial part of the magnetic field. 



3.4. The initial value and quasi-normal mode formulations 

spherical harmonics: In order to numerically integrate the system of equations we must split them using 
spherical harmonics, which removes the tensorial nature of the equations, and turns Eqs. (3-24) into algebraic 
relations. In the Appendix we have given an overview of the spherical harmonics we use, which were developed 
in Clarkson and Barrett (2003). 

A spherical harmonic decomposition of all variables then implies, from Eqs. (3-24), that the spherical harmonic 

jjtmUi Vsg'Vi V&V) i l »V> 
components of each of the variables iv , ^ a , M a , J a are proportional to the harmonic components 

of Xa ■ So, for example, for each I: 

*T = ^jxT\ (3-39a) 

f <4».<9> Lfijlg (L + Lg - L®) Jt&ig) 

(Lfij + L) Lg - (L - 



with similar relations for K 1 ^'^ , Mj^ 3 , J<f ■ 

Because the equations are linear in the second order variables, when we split into spherical harmonics, the 
equations decouple into two distinct subsets of opposing parity; the parity mixing which occurs between the 
magnetic field and the GW is contained in the interaction variables. We will call the set of equations containing Sy 
the even parity equations, and those containing §y the odd parity equations. Unfortunately, all the other variables 
are of the 'opposite' parity to E a in each system of equations, so this may cause confusion (so, e.g., fiy and xv are 
of even parity, etc.). 

INITIAL value FORMULATION: A useful form of writing ME is as wave equations. For this we will use the variables 

W s = r 2 Ss (even), (3-40a) 
W s = ^r 3 (3 s (odd). (3-40b) 

Then these variables satisfy the wave equations for each I: 

-W + W + AW-^W = S{ X a) (3-41) 



-13- 



where W — { W$ , W$ } and we have denned the even and odd source terms as 



OO OO 



5 S = -2Lr ^ ^ ^ Xlv + 2 X2v + 2 gV ' ; x 3v . (3-42a) 

(Leg — L) — (A® + L) L s I 



OO OO ( 



. 2 ..Ms*.*) l g (L ~ 3L,gg) ( % 



=1^ =2 



(£,^ - 4/ g - L) </> 2 r 2 + 4L (^) 2 £^ + Z g -L (w,) 1 
+ ~L~~L~$g X3 v +2r L-L* X4 v J' (3 " 42b) 

These wave equations may replace ME, and, more importantly, arc in the form of an initial value problem. We then 
have, for each parity, one forced wave equation for the EM field, plus two evolution equations for the interaction 
variables (xi and X3 ), plus two constraints (propagation equations); the set of four DEs for the interaction variables 
may be easily turned into a set of two coupled wave equations instead by eliminating two Xa variables (either xi 
and X3 or X 2 and Eliminating X2a using the first equation of (3-15) and X4a using the third equation turns 
the remaining two into wave equations: for the odd parity equations, we find 

-Xi v + Xi v +A Xlv =-2(<p + A)xi v +2x 3v 

+ [6£ \4> 2 -A 2 + (L„ L B ) r- 2 ] xT^ + W + 2A) X^ , (3-43a) 
~X 3v +X 3v +Ax 3v =1Lsgr Xi v -40X 3v 

+2L®r- 2 (cf, + A) xt^ + [7£ - 30 2 + {L g - L B ) r~ 2 ] ^ , (3-43b) 

with identical equations for the even variables. 

The full solution for the induced EM radiation is given by the variables W: for even perturbations, §y is given 
by (3-28a) and /3 V by (3-37); for odd perturbations, # v by (3-34) and (3y by (3-35). 



QUASI-NORMAL mode formulation using temporal harmonics: While the covariant equations above are 
given with time derivatives, allowing the problem to be put in the form suitable for solving as an initial value 
problem, it is often advantageous to use time harmonics. In particular, the effect of BH ringdown is conventionally 
studied by this method, as the ringdown phase is characterised by a set of quasi-normal frequencies (Nollert 1999; 
Kokkotas and Schmidt 1999), which are independent of the initial perturbation. We achieve this by replacing 
all dot-derivatives by a factor of iiv, with the usual understanding that subsequent equations are then for the 
spatial parts only (Clarkson and Barrett 2003), although formally it is significantly more complicated (Nollert 
1999; Kokkotas and Schmidt 1999; Andersson 1997). The harmonic function oj is defined with respect to the 
proper time, r, of observers travelling on u a , and satisfies Eq. (3-18); a is the constant harmonic index associated 
with time, t, measured by observers at infinity. Note that they are related by lot — at. 

The time derivative of the second order interaction variables J^ lm,lg) ^ \J> a ) ^ a ,M a ,J a acts only on 
the GW part of the term, because the time derivative of the magnetic field is already second order. Therefore, 
when these terms are split into time harmonics, and the interaction equations (3-23) are solved, the usual boundary 
conditions on the GW variable W a b will take effect - that the GW cannot propagate out of the horizon, or in from 
infinity. This implies that the allowed frequencies a must be discrete with positive imaginary part (Nollert 1999; 
Kokkotas and Schmidt 1999). This represents modes which decay exponentially in time, but whose amplitudes 
grow exponentially with radius. 

Our method presented here, which sets up the equations as a set of purely second order, linear, gauge-invariant 
differential equations, means that when we solve them wc don't view quadratic first order effects as quadratic forcing 
terms in the second order equations, but as second order quantities in their own right; the first order equations are 
forgotten about. Therefore the propagation equations governing the interaction variables, Eqs. (3-15) and (3-23), 
also must be confined to these frequencies. Hence, the coupling between the equations for the induced EM field 
and the interaction variables implies that the allowed independent frequencies of the induced EM radiation must 
be identical to those of the forcing GW - the quasi-normal frequencies; that is, the GW and EM radiation satisfy 
the same dispersion relation, and are in resonant interaction. Other frequencies correspond to EM waves which are 
not induced by the interaction terms with these boundary conditions (and form part of the homogeneous solution 
for the EM field); there is no need to consider these here. Therefore, when we split the system of equations using 
the time harmonics, each £ g picks out a set of allowed frequencies in the interaction equations, thus removing the 
summations over £ g in ME. For each £ g there is one system of equations for each quasi- normal frequency lo (1 ^ 

associated with that particular £ g . The complete solution for Sy, for example, may then be written schematically 
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for each I as 



<&=E«mo = E E 

£ g =2 <? 9 =2 U £{ U9 } n 



(3-44) 



where {^ 3 } n denotes the set of all quasi- normal frequencies for a given 

From the wave equations give above, it is clear that for each parity, while there are three EM variables, there 
are only two degrees of freedom in the EM radiation; in the even case for example these are Ws and Ws, resulting 
in a straightforward wave equation. We can of course stick to these variables in the QNM formulations of the 
problem, but the system is naturally first order in the variables Ssi and /3v (or d>s) in the even case once the extra 
degree of freedom is removed (similarly for the odd case) . There doesn't seem much advantage whichever way we 
choose the variables so we will remove the scalar (radial) parts of the EM field from the system of equations, using 
Eqs. (3-34) and (3-37). Our key equations then become: 



EVEN PARITY : 



\-(l + Lu 2 r 2 ) 0Xi v +2 X 



2v 



+ 4- 



Ig (L + Lg — Lc 



L) 



{Lf, 



L)L, 



-x 3v 



a n i r. a\ a In ^ \ n I 2 i -2\ J'*^ , J e &V n J e » 

/3v = - (|0 + 2A) A/ - wV v + E [ 2 ( w - l 9 r ) Xi v - ^X 3v - 2 X4v 



ODD PARITY : 



oo 



X 3v 



A/ = - (|0 + 2^) A/ - (-co 2 + Lr" 2 ) # v 
2 (L — 3L.og) l g 



+ E 



(L,j-L)r 2 



Xi v 



3v 



2X 



4v 



(3-45a) 

(3-45b) 
(3-45c) 

(3-45d) 

(3-45e) 
(3-45f) 



for each t g and u) G {^ g } n - Each parity consists of a set of six coupled ordinary differential equations in the radial 
parameter p. 



4. Numerical examples 

We have now set up the equations as a gauge-invariant linear system of differential equations in purely second 
order variables, in two different ways. The first is as a set of three coupled wave equations (for each parity), which 
may be numerically integrated as an initial value problem once some initial data is specified. The second is as a 
six-dimensional system of first-order ordinary differential equations (for each parity) which are Fourier decomposed 
in time, which is suitable for integration once appropriate boundary conditions are satisfied. There are of course 
advantages and disadvantages to both, which we discuss presently. 

While it would be desirable to be able to integrate these equations in a situation which is astrophysically 
accurate in some sense, this is quite a non-trivial problem as it involves specifying initial data from a fully non- 
linear integration of the field equations in a situation such as, for example, BH-BH merger. This is beyond the 
purpose of the present discussion, as we would like to get an overall estimate of the strength and importance of the 
effect in this first instance. 

In general, the summations over £ g and £ag in the equations for the generated EM radiation mean that these 
coupled systems of equations arc infinite dimensional. However, for a static magnetic field around a BH the 
dominant contribution to the field strength will be dipolar, and the GW emitted by a compact object will typically 
be dominated by the quadrupole radiation (for example, when two BHs collide head-on from an initially small 
separation, the emitted radiation is pure quadrupole (Price and Pullin 1994); other studies with high energy 
collisions support this conclusion (Cardoso et al. 2003; Cardoso and Lemos 2003a, b, 2002)). Therefore, in this 
section we will investigate numerically the £ g — 2, tag — \ interaction while ignoring the contribution from the 
others. 
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As we mentioned earlier, the case of an £gg = 1 magnetic field has two solutions, one which is uniform at 
infinity, and one which falls off at large distances like 1/r 3 , a dipole. Both of these are of interest astrophysically, 
as magnetic fields surrounding compact objects can extend considerable distances when supported by a plasma 
(i.e., a BH 'embedded' in an external magnetic field), but be purely dipolar close in. It is clearly important to 
distinguish the two cases in the Xa variables when we integrate the equations, in order to determine which type of 
field is responsible for what. For both solutions, the ratio SSy / ' 3§s is a known function of r, with no dependence 
on any boundary conditions - see, e.g., Eq. (3-6). This implies that 



Xi 

X3 



X2 
XA 



S8y 



uniform or dipole solution 



(4-1) 



where x% (* = 1> ' ' ' >4) represents either the odd or even parity part of Xia- The ratio SSy / 33$ is given by Eq. (3-6) 
for the dipolar field, while for the field which is uniform at infinity (characterised by SSc, — 0) it equals \<\>r. Thus, if 
we desire the magnetic field to be one of these solutions, we can use Eq. (4-1) to constrain the boundary conditions, 
or simply replace X3i Xi m t ne equations. In Fig. 1 we show a plot of the ratio — Xz/xi f° r the pure dipole field 
which shows how the dominant contribution to the interaction terms at large distances and close to the horizon is 
dominated by xi (or Xn\ the figure for — x^lXi is identical). Thus, X3 an d X4; containing the radial part of the 
magnetic field, only contribute significantly in the vicinity of the photon sphere. We will consider only the pure 
dipole solution here, and hereafter remove X3,X4 using Eq. (4-1) (we remove these two because, as Fig. 1 shows, 
replacing xi ; Xi would make numerical solutions become unstable at small and large distances) . 
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Fig. 1. — For a pure dipole background magnetic field this figure, which is a plot of Eq. (4-1), shows how the 
relative contributions from the interaction terms are dominated by \i (or X2)> except in the region just inside the 
photon sphere (the peak is at about r ~ 2 Am), where X3 becomes significant (recall that X3 is defined using the 
angular gradients of the radial part of the background magnetic field). The 'tortoise' coordinate r* is defined by 
Eq. (4-2) below. 



The induced EM radiation will of course be of much higher amplitude far from the BH if we allow for the 
presence of the uniform magnetic field as part of the static background, as the interaction distance will be vastly 
increased. For the pure dipole magnetic field, the interaction distance is effectively curtailed at large r, because 
the magnetic field strength falls off so fast that S ^> x by r ~ 20m or so. In astrophysical situations where the 
magnetic field extends far from the source (supported by an accretion disk, or entangled in the ejected envelope of 
the progenitor star, for example), we would expect further, linearly growing amplification, (beyond say r ~ 20m) 
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over the amplification we report below. This will be studied at a later date when a plasma is included into the 
discussion, but should be borne in mind in what follows. 

Hereafter we shall set m = 1 (which just defines the units of r), and we shall use the tortoise coordinate r* of 
Regge and Wheeler (Regge and Wheeler 1957) defined by 



Because the system of equations we are investigating are linear, the units we use are physically irrelevant, and is 
tied into the physical amplitude of our initial data which we normalise at unity (so that if units are chosen for % 
say we can immediately read off the actual amplitudes for $). 



Here we envisage the following situation: at some initial time t = the interaction is 'turned on' with some 
typical initial profile for the GW [i.e., the tensor W a b, which translates in this case to X a (t = 0) = Xo]i & t which 
time the induced EM field is zero, but with non-zero second time derivatives ('acceleration'). Although intuitively 
reasonable for modeling a situation such as BH formation or where the magnetic field becomes very strong very 
quickly, say, we require this switching on of the interaction because otherwise the ME will not be consistent for a 
general Xo- 

A common way of specifying initial data for this type of problem is to consider GW scattering off a BH, with 
the initial data given by a static narrow Gaussian peak at some distance from the hole (Andersson 1997). This 
then splits in two as the RW equation is evolved, with the part falling into the hole of most interest: this scatters 
off the photon sphere and starts the black hole vibrating (roughly speaking), with a characteristic waveform 
which is largely independent of the initial data, dominated by the quasi-normal modes of the BH (which only 
depend on its mass) (Andersson 1995, 1997; Sun and Price 1998; Nollcrt 1999). We will use this scenario with 
Wo ~ cxp (— (r* — 20) 2 ) at t = 0, which we normalise so that at t = and r* = 20, xi v = 1- We will not consider a 
pulse originating further from the hole because the dipole field falls off so fast with distance; the qualitative results 
remain the same. 

We then evolve our key equation (3-41) and the wave equation for xi [modified by replacing X3 with Eq. (4-1) 
as discussed above] with this initial data. This then gives the solution for Ws, which we convert to Sy. Results 
are shown in Figs. (2) and (3) for log 10 \S\j\- These figures show the EM radiation generated and subsequently 
amplified during the scattering of the GW off the photon sphere. The ringing of the BH then generates a continuous 
stream of EM radiation, which at its peak is over two orders of magnitude larger than the initial pulse of radiation 
(by the time it is reflected back out to r* = 20). This radiation mirrors very closely the GW waveform making it 
a suitable EM counterpart for GW emission. 



We shall now integrate the equations in the frequency domain, summing over the QNMs of the BH, which will 
tell us about the strength of the interaction in the latter stages of a perturbation of a BH independently of the initial 
perturbation (Andersson 1997). We imagine that the interaction starts at t = at some inner radius ro, so for r < r 
we assume that £$ = /?s = 0, while x does its own thing; at r — r we choose boundary conditions for each w n such 
that all EM terms and their derivatives are equal to zero; for want of accurate boundary conditions for the GW, we 
randomly 4 choose x^ n>> — Xo- In order to compare differing amplifications for each parity, we use the same Xo for 
both parities. We then integrate Eqs. (3-45a) - (3-45) out to some r — r max for each QNM frequency ui n . Then, 
for each variable at r — r m ax we can simply add up the QNMs. This then gives a good approximation to the time 
decay of the signal as it passes r = r m &x after t > tmax = ^max — f o + 2m In [(r m ax — 2m) / (ro — 2m)] (Andersson 
1997; Nollert 1999). We use the first twelve QNM frequencies as tabulated in Nollcrt and Schmidt (1992) for 
o~ n = \<j)ru n - see Eq. (3-18). 

In Figure 4 we show a typical result of this integration, for an observer situated at r ~ 65, with ro = 2.05. 
The generated electric field is shown, for both parities, as is the largest interaction variable, which is xi v m this 
case. The units of the graph are arbitrary: dividing each variable by |xo| (to make each variable dimensionless) , 
say, will merely shift all the curves up or down. At large distances from the source, the behaviour of the fields can 
be represented as an amplitude over a potential of the distance function. In the case of the gravitational wave, 
the fall-off scales like 1/r, while for a spherical electromagnetic wave it behaves as 1/r. At the same time, the 
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4.1. The initial value problem 



4.2. The quasi-normal mode approach 



4 Although this may seem somewhat arbitrary, it is no more arbitrary than choosing a Gaussian distribution as in the last section. 
We have performed the numerical integration below for many different choices of xo> an d the results are qualitatively similar. 
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Fig. 2. — The induced EM radiation from static Gaussian initial data for the GW at t = 0, such that the EM 
field is zero. This pulse splits into two, one falling into the hole and the other propagating to infinity. The part 
falling into the hole is partially reflected at about t ~ 15 generating the 'ringing' we see later at t = 50 in the thick 
curve (modulated by the magnetic field here - this is Xi v )- During this in- fall, the GW-B interaction produces 
substantial amounts of EM radiation which is reflected back away from the hole, and is further increased by the 
subsequent BH ringing. 

background magnetic dipole field has l/r 3 -dependence. We can therefore normalise with the respect to the fall-off 
of the field strengths, in order to get a scale-invariant form of the amplification. In the situation given in Fig. 4, 
normalising the curve for xi v raises it up by 31og 10 (65/2.05) ~ 4.5. Hence, at large distances from the source the 
scale-invariant amplification of the EM radiation is over two orders of magnitude larger the the magnitude of the 
GW times the magnetic field strength. At this distance from the source the interaction is no longer taking place 
to a significant degree, implying that this level of amplification is a generic feature. Note also from Fig. 4 that the 
amplification of the electric field parities is different, implying that the EM wave is polarised. 



To estimate the implications of this amplification we have found, consider the case of a compact object such as 
a BH or neutron star. The interaction between gravity waves and the magnetic field is quantified by the variable 
Xi ~ hiBi, where hi is the amplitude of the gravitational wave at the onset of the interaction, Bi ~ B s (r s /ri) 3 
is the field strength at the distance r,; from the compact object with 'radius' (e.g., surface of a neutron star or 
BH horizon) r s and surface magnetic field B s . The interaction produces an EM signal E ont which is typically two 
orders of magnitude larger at r cut _ ff than the original perturbation \i a t ^cut-off) where the interaction effectively 
switches off. From Fig. 3 we extrapolate that at r cut _ ff ~ 40m we have roughly Xcut-off ~ 10 -2 ' 5 Xi? leading to an 
induced electric field strength 



4.3. Estimates 




(4-3) 
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Fig. 3. — The induced EM radiation from static Gaussian initial data for the GW at t = 0, such that W ~ 
exp (— (r* — 20) 2 ), normalised so that at t = and r* — 20, xi v = 1; at t = the EM field is zero. This pulse splits 
into two, one falling into the hole and the other propagating to infinity. The part falling into the hole is partially 
reflected generating the 'ringing' we see at t = 55 (modulated by the magnetic field here in xi v ), in the wave which 
is moving to the right . During this in-fall, the GW-B interaction produces substantial amounts of EM radiation, 
which is further increased by the subsequent BH ringing. Thus we see at t = 55 the EM field generated is two 
orders of magnitude larger than the initial pulse, and over three orders of magnitude larger than the interaction 
terms at the same distance from the hole. At t = 80 we see that the main amplification has taken place with the 
peak of the induced EM radiation moving past r* ~ 40. 



The induced signal attenuates inversely with distance D outside the interaction region, e.g. r > 40m. At a distance 
D from the source, the spectral energy flux $ w = (^eocE 2 )/uj can then be calculated to be 

As an example, consider a magnetar of mass m — 1.5 Mq with radius r s = 9 km, e.g. twice its Schwarzschild 
radius, and take r-j = 4 r s and r cut-off = 90 km, assuming a magnetic field strength B s in the range of 10 5 to 10 10 
T. An occurring instability such as a supernova explosion or a bar mode instability is likely to produce a GW with 
hi ~ 10~ 3 and frequency uj of about 1 — 10 kHz (Andersson 2003), which is also the frequency of the induced EM 
wave. This leads to E ou t ~ 5 x 10 8 — 10 13 Vm _1 . If such an event happened within our galaxy (D ~ 10 kpc), 
<f> w ~ 10 5 — 10 15 Jy, if we assume that 10% of the signal's energy undergoes mode conversion shifting the frequency 
up to 30 — 300 kHz (Marklund, Brodin, and Dunsby 2000). To achieve a higher detection rate, one has to gather 
events from a farther distance. Events within the Virgo Cluster (D ~ 15 Mpc) would have flux $/ ~ 10 _1 — 10 9 
Jy. The proposed radio telescope Astronomical Low Frequency Array (Jones et al. 1998) is expected to operate 
in the range from 30 kHz to 30 MHz, with minimum detection level of 1000 Jy, making such events an exciting 
possibility for indirect gravitational wave detection. 



5. Conclusions 

We have investigated the scenario of GW around a Schwarzschild BH interacting with a strong, static, magnetic 
field. This interaction produces a stream of EM radiation mirroring the BH ringdown, with a stronger amplitude 
than one may expect from estimates of the interaction in flat space, due to non-linear amplification in the vicinity of 
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Fig. 4. — Late time behavior showing the temporal evolution of the generated EM waveforms for an observer at 
r = 65771, using a quasi-normal mode expansion. Choosing boundary conditions at r = 2.05m, the QNM sum 
will approximate the true solution for this observer after t ~ 76m. We see that the EM waveform is substantially 
larger than the interaction variables, for which we have shown the largest, xi v - The differing amplification of the 
electric field for each parity is also not insubstantial, with the even parity case being larger than the odd, in this 
case, demonstrating polarisation of the induced EM radiation. 

the photon sphere. This interaction may play an important role in GRBs and perhaps some SN events, in addition 
to neutron star physics, and may be a useful mechanism to aid in GW detection. 

We converted the Einstein-Maxwell equations into a linear, gauge-invariant system of differential equations by 
utilising the 1+1+2 covariant approach to perturbations of Schwarzschild. We also introduced a set of second order 
'interaction' variables to aid in simplifying the derivation, and a new variable for the magnetic field, both of which 
made the system of equations manifestly gauge-invariant. It was then a simple matter to convert the system of 
equations into wave equations for integration as an initial value problem, or as a harmonically decomposed (in time) 
system of first-order ordinary differential equations, which could then be integrated using a BH quasi-normal mode 
expansion, an important approximation method for late time behaviour. We integrated the system of equations 
using both of these techniques. 

A key point of this paper was to set up a suitable formalism to study this GW-B interaction around a BH, and 
to put the equations into a suitable gauge-invariant form for numerical integration. The next step is to include a 
plasma, as various plasma instabilities could be induced by such process, making detection of this sort of induced 
radiation a genuine possibility. This will also help model some of the relativistic effects which take place after a 
SN explosion. In fact, EM waves in a plasma in an exact Schwarzschild spacetime are pretty complicated and 
unexpected (Daniel and Tajima 1997), so it is an interesting question in its own right to ask what happens when 
GW are thrown into the mix. 

This project was financially supported by NRF (South Africa) and Sida (Sweden). CAC would like to thank 
the Department of Electromagnetics at Chalmers University of Technology for their hospitality while part of this 
work was undertaken. We would like to thank Richard Barrett and Marco Bruni for useful discussions. 

A. Spherical harmonics 

Here we briefly review the spherical harmonic expansion, developed in Clarkson and Barrett (2003) appropriate 
to the formalism, for easy reference. These allow us to remove all (^-derivatives from the equations. Note that all 
functions and relations below are defined in the background only; we only expand second-order variables, or first- 
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order variables which form part of a quadratic second-order variable so zeroth-order equations are sufficient. 

We introduce spherical harmonic functions Q — Q( e > m \ with m = —£, ■ ■ ■ ,£, defined on the background, such 

that 

S 2 Q = -£(e+l)r- 2 Q, Q = = Q. (Al) 

We also need to expand vectors and tensors in spherical harmonics. We therefore define the even (electric) parity 
vector spherical harmonics for £ > 1 as 

QP=r5 a QW => Q a = = Q a , 5 2 Q a = (1 - £ {£ + 1)) r - 2 Q a ; (A2a) 

where the (£) superscript is implicit, and we define odd (magnetic) parity vector spherical harmonics as 

Q^=re ab 6 b QW Q a = = Q a , 8 2 Q a = (I -£{£+!)) r - 2 Q a . (A2b) 

Note that Q a — e a bQ b <S=> Q a = —e ab Q b , so that e ab is a parity operator. The crucial difference between these two 
types of vector spherical harmonics is that Q a is solenoidal, so 

5 a Q a = 0, while S a Q a = -£(£+l)r- 1 Q. (A3) 

Note also that 

e ab S a Q b = 0, and e ab 6 a Q b = £ (£ + 1) r^Q. (A4) 

The harmonics are orthogonal: Q a Q a — (for each £). Similarly we define even and odd tensor spherical harmonics 
for £ > 2 as 

Q ab = r 2 5 {a S b} Q, => Q ab = Q = Q ab , S 2 Q ab = [<j> 2 - 3£ - I (£ + 1) r" 2 ] Q ab ; (A5a) 

and 

Qab = r 2 e c{a S c 6 b} Q, => Q ab = = Q ab , 8 2 Q ab = [0 2 - 3£ - £ (£ + 1) r~ 2 ] Q ab , (A5b) 

which are orthogonal: Q ab Q ab — 0, and are parity inversions of one another: Q ab — —£ c {aQ b } ^ Qa b — £ c {aQ b f- 
We can now expand any second-order scalar * in terms of these functions as 

oo m— £ 

* = E ^f m) Q {t ' m) - *sQ, (A6) 
e=o m=-e 

where the sum over £ and m is implicit in the last equality. The S subscript reminds us that ^ is a scalar, and 
that a spherical harmonic expansion has been made. Due to the spherical symmetry of the background, m never 
appears in any equation so we can just ignore it. Any second-order vector ^ a can now be written 

oo 

*a=E + *V = *vQa + * V Qa. (A7) 

Again, we implicitly assume a sum over £ in the last equality, and the V reminds us that \l/ a is a vector expanded 
in spherical harmonics. Any second-order tensor may be also be expanded 

oo 

*ab = J2 ^Qab + V^Qab = *tQo4 + ^jQab- (A8) 

1=2 

Further useful identities are to be found in Clarkson and Barrett (2003). 
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